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Abstract 

The non-linear evolution of the halo population is followed by solving the continuity equation under 
the hypothesis that haloes move by the action of gravity An exact and general formula for the 
Eulerian bias field of dark matter haloes in terms of the Lagrangian bias is expanded at second- 
order including the presence of primordial non-Gaussianity Particular attention is paid in defining a 
gauge-invariant bias which is necessary when dealing with relativistic effects and measured quantities. 
We show that scale-dependent effects in the Eulerian bias arise both at first- and second-order 
independently from the presence of some primordial non-Gaussianity. Furthermore, the Eulerian 
bias inherits from the primordial non-Gaussianity not only a scale-dependence, but also a modulation 
with the angle of observation when sources with different biases are correlated. 

1 Introduction 

Cosmological inflation [1] has become the dominant paradigm to understand the initial conditions for 
the Cosmic Microwave Background (CMB) anisotropies and Large Scale Structure (LSS) formation. 
This picture has recently received further spectacular confirmation by the Wilkinson Microwave 
Anisotropy Probe (WMAP) seven year set of data [2]. Present and future [3] data may be sensitive 
to the non-linearities of the cosmological perturbations at the level of second- or higher-order per- 
turbation theory. The detection of these non-linearities through the non-Gaussianity (NG) [4] has 
become one of the primary experimental targets. 

A possible source of NG could be primordial in origin, being specific to a particular mechanism 
for the generation of the cosmological perturbations. This is what makes a positive detection of 
NG so relevant: it might help in discriminating among competing scenarios which otherwise might 
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be indistinguishable. Indeed, various models of inflation, firmly rooted in modern particle physics 
theory, predict a significant amount of primordial NG generated either during or immediately after 
inflation when the comoving curvature perturbation becomes constant on super- horizon scales [4]. 
While standard single-field [5] and two(multi)-field [6] models of inflation generically predict a tiny 
level of NG, 'curvaton-type models' [7, 8, 9], in which a significant contribution to the curvature 
perturbation is generated after the end of slow-roll inflation by the perturbation in a field which has a 
negligible effect on inflation, may predict a high level of NG [10]. Alternatives to the curvaton model 
are those models characterized by the curvature perturbation being generated by an inhomogeneity 
in the decay rate [11, 12] of the inflaton field. Other opportunities for generating the curvature 
perturbation occur at the end of inflation [13] and during preheating [14]. All these models generate 
a level of NG which is local, since the NG part of the primordial curvature perturbation is a local 
function of the Gaussian part generated on superhorizon scales. It has now become common to 
parametrize the level of NG through a dimensionless quantity /nl which sets the magnitude of the 
three-point correlation function [4]. In momentum space, the three point function (bispectrum), 
arising from the local NG is dominated by the so-called "squeezed" configuration, where one of the 
momenta is much smaller than the other two and it is parametrized by the non-linearity parameter 
/^l- Other models, such as DBI inflation [15] and ghost inflation [16], predict a different kind of 
primordial NG, called "equilateral", because the three-point function for this kind of NG is peaked 
on equilateral configurations, in which the lengths of the three wave- vectors forming a triangle in 
Fourier space are equal [17]. The equilateral NG is parametrized by an amplitude f^ 11 [18]. Present 
limits on NG are summarized by -10 < f§£ < 74 and -214 < /^ uil < 266 at 95% CL [2, 19, 20]. 

It is clear that detecting a significant amount of NG and its shape either from the CMB or from 
the LSS offers the possibility of opening a window into the dynamics of the universe during the 
very first stages of its evolution and to understand what mechanism gave rise to the cosmological 
perturbations. Besides in the CMB anisotropics, NG is particularly relevant in the high-mass end 
of of density perturbations, i.e. on the scale of galaxy clusters, since the effect of NG fluctuations 
becomes especially visible on the tail of the probability distribution function [21]. Furthermore, 
and more relevantly for us, primordial NG also alters the clustering of dark matter halos inducing 
a scale-dependent bias on large scales [22]. According to the peak-background split theory [23] 
the underlying idea behind the generation of a local bias is that galaxies tend to form in regions 
where the dark matter density field is larger than some threshold value in Lagrangian space. The 
collapse of objects on small scales is ascribed to the high frequency modes of the density fields, while 
the action of large-scale structures of these non-linear condensations is due to a shift of the local 
background density. As primordial NG generates a cross-talk between short and long wavelengths, 
it alters significantly the local bias and introduces a strong scale dependence in it. As a result, 
measuring the clustering properties of haloes is a sensitive probe of primordial NG which could be 
detected or significantly constrained by the various planned large-scale galaxy surveys, both ground 
based (such as DES, PanSTARRS and LSST) and in space (such as EUCLID and ADEPT) [24]. 

When analyzing the impact of NG onto the bias of dark matter haloes various points should be 
addressed. Being the effect computed at second-order in perturbation theory one should consistently 
calculate how to go from the Lagrangian bias to the Eulerian bias at the same order in perturbation 
theory. Furthermore, since the primordial NG manifests itself on large cosmological scales one should 
treat carefully the relativistic effects. This automatically calls for a gauge-invariant formulation of the 
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observables at hand. In this paper we address these points and show that a refined gauge-invariant 
treatment of the Eulerian bias at second-order in perturbation theory leads to the prediction that 
the bias is scale-dependent on large scales even in the absence of primordial NG and that the latter 
generates an angular modulation if sources with different biases are cross-correlated. 

The paper is organized as follows. In section 2 we describe how to obtain the Eulerian bias 
description in terms of the Lagrangian bias description at the non-perturbative level. In section 
3 we specificaly deal with the gauge-invariant formulation and in section 4 we proceed with the 
computation of the gauge-invariant Eulerian bias at second-order in perturbation theory. Finally, 
section 5 contains our conclusions. 

2 Eulerian description of the bias from the local La- 
grangian bias 

In this section we describe how to obtain an expression for the Eulerian halo-to-mass bias starting 
from the local Lagrangian bias. According to the local Lagrangian description, the sites of the galaxy 
formation are identified with specific regions of the primordial density field. The primordial galaxy 
density field measuring the (smoothed) overdensity of galaxies in fieri at the Lagrangian position q 
at a given time r <C 1 is biased with respect to the primordial (linear) matter CDM density at the 
same location and at the same time. The Eulerian bias is obtained by integrating the continuity 
equations for the mass and for the halo number density where galaxies are supposed to reside. The 
procedure follows and generalizes the one provided in Rcf. [25] for the Newtonian case. We consider 
a spatially flat Universe filled with a cosmological constant A and a non-relativistic pressureless 
fluid of Cold Dark Matter (CDM), whose energy-momentum tensor reads T^ v = pu^Uy where 
{u^u^ 1 = —1) is the comoving four- velocity. 

Following the notations of Ref. [26], the perturbed line element around a spatially flat FRW 
background reads 

ds 2 = a 2 (r){-(l + 2<j))dT 2 + 2u i dTdx i + [(1 - 2V>)<% + XijWdx 1 } . (1) 

where a(r) is the scale factor as a function of conformal time r. Here each perturbation quantity 
can be expanded into a first-order (linear) part and a second-order contribution, as for example, the 
gravitational potential 4> = <fi^ + 4>^ /2. Up to now we have not choosen any particular gauge. We 
can employ the standard split of the perturbations into the so-called scalar, vector and tensor parts, 
according to their transformation properties with respect to the 3-dimensional space with metric 
6ij, where scalar parts are related to a scalar potential, vector parts to transverse (divergence- free) 
vectors and tensor parts to transverse trace-free tensors. Thus <j) and ip, the gravitational potentials, 
are scalar perturbations, and for instance, Cj^ = diU^ where oj^ is the scalar part and oj\ 

is a transverse vector, i.e. d l uj^ = ((r) = (1,2) stand for the rth-order of the perturbations). 
The symmetric traceless tensor Xij generally contains a scalar, a vector and a tensor contribution, 
namely Xij = ^ijX + diXj + djXi + Xij> where Dij = didj — (l/3)V 2 <5jj, Xi ls a solenoidal vector 
{d l Xi = 0) and Xij represents a traceless and transverse (i.e. d l Xij = 0) tensor mode 1 . The spatial 

1 ln what follows, for our purposes we will neglect linear vector modes since they are not produced in 
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projection tensor orthogonal to the fluid velocity u M is defined by 

V> = 9tiv + u^uy, {KK = K u » = °)- ( 2 ) 

It is also useful to introduce the familiar decomposition 

V v u^ = o^ v + to^u + -Qh^ - a^Uy, (3) 

where we have defined the (symmetric) shear tensor a^ u , the (antisymmetric) vorticity tensor u^, 
the volume expansion scalar 6 = V M ii M and the acceleration a M = u u V u u^. Notice that reduces 
to 37i (H being the Hubble rate in conformal time) in the homogeneous and isotropic FRW case. 
Our starting point is the conservation of the energy- momentum tensor of the CDM fluid, 

= 0, (4) 

which yields, after projecting along u v , the continuity equation valid at any order in perturbation 
theory 

p(x,r) + e(x,r)p(x,r) = 0, (5) 

where the dot indicates differentiation along u^, that is p = u^V^p. If we now assume that a halo 
population of mass M and formation time Tf is conserved in time and evolves exclusively under the 
influence of gravity with an unbiased velocity 2 , meaning that the CDM fluid and haloes are moving 
with the same velocity, its number density ph( x > T ) = Ph( x > t\M, Tf ) has to satisfy the continuity 
equation 



p h (x,T) + 6p h (x,r) = 0. (6) 

Notice that Eqs. (5) and (6) are non-perturbative and valid in any gauge. This is welcome as we want 
to derive the Eulerian bias factor which, being a physical observable, must be a gauge-independent 
quantity. While comparing the theoretical predictions (the matter power spectrum) with obervations 
(the galaxy power spectrum) does not represent a problem on sub-horizon scales where the matter 
density perturbations computed in the different gauges all coincide, it is a delicate operation on 
scales comparable with the horizon where different gauges provide different results even at the linear 
level (see, e.g., [28]). Truly gauge-independent perturbations must be exactly constant in the back- 
ground spacetime. This apparently limits ones ability to make a gauge-invariant study of quantities 
that evolve in the background spacetime, e.g. density perturbations in an expanding cosmology. In 
practice one can construct gauge-invariant definitions of unambiguous, that is physically defined, 
perturbations (see, e.g., the discussion of Ref. [29]). These are not unique gauge-independent per- 
turbations, but are gauge-invariant in the sense commonly used by cosmologists to define a physical 
perturbation. There is a distinction between quantities that are automatically gauge-independent, 
i.e., those that have no gauge dependence (such as perturbations about a constant scalar field), and 
quantities that are in general gauge-dependent (such as the curvature perturbation) but can have a 
gauge-invariant definition once their gauge-dependence is fixed (such as the curvature perturbation 

standard mechanisms for the generation of cosmological perturbations (as inflation), and we also neglect 
tensor modes at linear order, since they give a negligible contribution to LSS formation. 

2 This approximation is accurate if one is interested, as we are, on the bias at large-scales [27]. 
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on uniform-density hypersurfaces). In other words, one can define gauge-invariant quantities which 
are simply a coordinate independent definition of the perturbations in the given gauge. This can be 
often achieved by defining unambiguously a specific slicing into spatial hypersurfaces. In this sense 
it should be clear that one may define an infinite number of, e.g., gauge-invariant density contrasts. 
Which one to use is a matter that can be decided only considering how the determination of a given 
observable is performed. We will come back to this point later. For the time being it suffices to 
say that, when expanded at a given order in perturbation theory, Eqs. (5) and (6) may be used 
to find the evolution of the gauge-independent (in the sense just described) CDM and halo density 
contrasts, 5 S1 and 5^. These quantities evolve following the same dynamics. This means that their 
density contrasts 5 S1 and at a given time will be related to their values at some initial time T m <C 1 
through the relation 

1+jgW) = l + ^ gi (x,r) 
l + «5f(q) l + ^ gi (q) ' 

In Eq. (7) by 5f*(q) = 5^(q\M, Tf) we mean the initial (Lagrangian) halo density fluctuation at 
some primordial time when x(r = r; n ) = q. Expanding Eq. (7) up to second-order and setting 

*h = 4 1) + 5 5 h ) . wefind 

5f W (x, r) ~ % {1) (q) + (>« ( X , r) - 5*™ (q)) , (8) 

and 

\€ m (x, r) c l -5f^ (q) + i (>(2) (x> T) _ 6 ,m (q) ) 

+ (q) - S*W (q)) (x, r) - (q)) , 

(9) 

The expressions (8) and (9) are the key relations which permit to relate the Lagrangian bias to 
the Eulerian one. Before continuing though, we come back to the issue of which gauge-invariant 
contrasts we should take. 



3 On the gauge- invariant formulation 

As we pointed out before, there is an infinite number of ways to define gauge-invariant density 
contrasts which differ by other gauge-invariant combinations. Since one observes galaxies rather 
than the underlying matter distribution and the latter at the source galaxy position is related to 
the mean matter density at the observed redshift z, a good choice to define gauge- invariant density 
constrasts related to each other by a bias factor seems the one involving the observed redshift z [28] . 
At first order a coordinate transformation reads x M — >• x^ — £^ where £^ = (af 1 ',^ 1 ' 1 ). The matter 
density contrast transforms as 5^ — > 5^ + p/pa^\ where now dot stands for differentation with 
respect to the conformal time and p ~ (1 + z) 3 is the background matter energy density; similarly 
the first-order perturbation of the observed redshift transforms as z^ —> z^ + zoS^ (here z is 
the unperturbed redshift). Going to the uniform redshift gauge where the linear perturbation of 
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redshift vanishes relates to the linear perturbation of redshift in the old gauge, = —z^/z. 
Therefore the gauge invariant definition of the matter density contrast (and similarly for the halo 
one) is [28] 

JS i(D = ^(1) _ 3^1 (10) 
l + z y ' 

At second-order the procedure is more involved, but straightforward. The coordinate transforma- 
tion reads x^ -> + £^ + \ (tfo^fa + £f 2 )) where £(2) = (« (2) , £ (2)i )- Under this coordinate 
transformation the density matter contrast and the redshift perturbation transform as 

P \P P P J 1 ' \P P J 

z m z( 2 ) + ia( 2 ) + a( 1 )(za( 1 )+id( 1 )+2i( 1 ))+^ 1) (i^a« + 2a^( 1 )). (11) 

Going to the uniform redshift gauge where second-order perturbation of the redshift vanishes gives 

_(2) r (iU(i) £L 

a m = _ z lL + z l^l_^ diZ w. ( i2) 

z z A z 

To completely solve the uniform redshift gauge at second-order we must also specify the first-order 
spatial gauge shift fl^ • A natural choice is to pick worldlines comoving with the fluid. The (scalar) 

velocity transforms as v^ 1 — > — £^ . Thus from an arbitrary spatial gauge we can transform to 
the comoving gauge by the spatial gauge transformation ^ = f vf^dr. Such a choice leads to the 
second-order gauge invariant matter density constrast (and similarly for the halo one) 



- 3^+(^-3^)/^> 

The next step amounts to determing the expression of the redshift perturbation in terms of the 
perturbations of the metric (1) and other quantities. Photons suffer a redshift z during their travel 
from the emitter 8 to the observer O; the emitted frequency uj£ and the observed one uiq are related 
by ojo = ws/{\ + z). Here to = —g^ v u^k v , where is the four-velocity of the observer or emitter and 
k u = dx v /d\ is the wave vector of the photon in the conformal metric, tangent to the null geodesic 
x v {\) (A is the affine parameter) followed by the photon from the emission to the observation point. 
We do not report the full calculation of the redshift perturbation which basically amounts to solving 
for the photon trajectory. The computation can be found in Ref. [30]. Expanding the frequency as 
uj = uj(1 + a/ 1 ) + ^uj^), at first-order one obtains 
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l + z 



(1) 

Xq 



,(!) 



-^ +v ^e t -v^e l + I 1 (X £ ), 



(14) 



where e* indicates the zero-th order three-dimensional vector indicating the photon direction from 
which they arrive at the observer O. In the expression above one recognizes the Sachs- Wolfe effect 
due to the change in the gravitational potential at the source' and observer's points, the Doppler 
contribution due to the peculiar velocities of the emitter and the observer and the integrated Sachs- 
Wolfe effect along the photon trajectory. 

At second-order expressions are more involved and the perturbation of the redshift may be 
written as [30] 



1 z 



(2) 



21 + z 



where 



+ 



1 

2 

2 <Po 



W ? ) - W g ) -2 W « W g> + 2(a;g>) : 



,(2)< 



(2)< 



2*g> 



_ T (i)o i(i) _ 



+ x 



(1)0, 



- i(*g')' 

v$ l ei - < 



^ +v£ )i e i + I 1 (\ £ ) 



(i)< 



1 



1 (l) (i)i . (i)i 



- 4V (-4 1)l + *§* + ^gV - xg } % - Il(Xs)) 
+ (fpei " "S** + A(A f ))* + 2 (*g> - ( 



(i)i 
w e i 



+ A(A £ )) 



^i(Ae) 
/ 2 (A f ) 



A £ 



/ ' dA -i( 2 ) - (4 W - xS^KifeW* + e l fc (1)0 ) + 2fc( 1 )°i( 1 ) 



A m = v ,( 2 ) + (2) + ^) e ,_i x g) e v 



(15) 



(16) 
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and 

k^°(X £ ) 
k^(X £ ) 



(l)i 



(i)i 



2V. 



UJ 



(Af-Ao) 2^ ) e' + cj^ / ' -xo 



+ 



dA 



ej-iiCAf), 



UK?, 



- ! " dX \2^e l + co Wi - x (1)ij e j + (A £ - A)^ 1 ^ 



(17) 



These long expressions allow the determination of the gauge-independent density contrasts. In 
practice, one evaluates them adopting the most appropriate gauge. One possible convenient choice 
for the determination of the Eulerian bias from the local Lagrangian bias is represented by the 
comoving-orthogonal gauge as we will explain in the following. Let us stress that the gauge-invariant 
expressions found in this subsection would be also useful if, instead of computing the Eulerian bias 
derived from the local Lagrangian bias, one adopts the bias model described in Ref. [31]. In this 
approach the (smoothed) galaxy number density field at a given position x and time r is assumed 
to be a local function of the (smoothed) underlying CDM mass density at the same location and 
instant, 5ff = bf5 gl + (<5 g1 ) 2 + • • • . This approach is essentially phenomenological and it is a 
priori devoid of any insight about the dynamics of the clustering. 



4 The computation of the gauge- invariant Eulerian bias 

In this section we proceed with the computation of the Eulerian bias adopting the gauge-invariant 
formulation of the matter density perturbation described previously. It is convenient to perform such 
a computation collapsing to the comoving-orthogonal gauge. Indeed, as we restrict ourselves to the 
case of irrotational dust plus A, the comoving-orthogonal gauge is also synchronous. Indeed, the fluid 
four- velocity field can be written as = (1/a, 0, 0, 0), so that q represents the comoving Lagrangian 
coordinate for the fluid element. The possibility of making the synchronous, time-orthogonal gauge 
choice and comoving gauge choice simultaneously is a peculiarity of fluids with vanishing spatial 
pressure gradients, i.e. vanishing acceleration, which holds at any time, i.e. also beyond the linear 
regime. 

The choice of the comoving-orthogonal gauge to evaluate the gauge-invariant Eulerian bias is 
motivated by various reasons. A simple analytic model for the gravitational clustering of dark 
matter haloes to understand how their spatial distribution is biased relative to that of the mass was 
developed in Ref. [32]. The statistical distribution of dark haloes within the initial density field 
(assumed Gaussian) is determined by an extension of the Press-Schechter formalism and is done 
therefore at the Lagrangian level. One then expects that the gauge-invariant Eulerian description 
to be therefore simpler to formulate through the synchronous gauge. Furthermore, the non-linear 
spherical collapse description necessary to compute the halo mass function through the (extended) 
Press-Schechter approach, requires the choice of the comoving-orthogonal gauge. 

In the synchronous and comoving gauge the line element reads 

ds 2 = a 2 (r) [-d 2 r + hij(*,T)djdx>] . (18) 
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Working in the synchronous comoving gauge, the spatial coordinate does not evolve with time, 
x(t) = x(r = r; n ) = q. From Eqs. (8) and (9) evaluated now in the comoving-orthogonal gauge, 
we may deduce the gauge-invariant Eulerian bias parameters. To simplify our expressions we limit 
ourselves to the case of a pure CDM-dominated universe. One can easily check that the gauge- 
invariant matter density contrast (10) expressed in the comoving-orthogonal gauge is given by 3 

<5Si(1)(q ' T) = 37^ V ^ (q) + ^ (q) + re ^( q )' ( 19 ) 

where <£>(q) is the peculiar scalar gravitational potential and we have got rid of the terms defined at 
the observer's point which may be absorbed in the monopole term. Notice that in Fourier space the 
time evolution of the gauge-invariant density contrast as 



^(D( k ,r) = / gi (k,r)«5W(k,r), 

/ g i(k,r) = 1-^-6^, (20) 

where 5^ (k, r) = — ^-£-ip{\s) is the matter density contrast in the synchronous gauge. The key 
point is now to calculate the halo field <5 gl (q). The approach to the clustering evolution is based on 
a generalization of the so-called peak-background split, first proposed in Rcf. [23] which basically 
consists in splitting the mass perturbations in a fine-grained (peak) component 5^ k filtered on a scale 
R and a coarse-grained (background) component 5® filtered on a scale i?o S> R. The underlying 
idea is to ascribe the collapse of objects on small scales to the high frequency modes of the density 
fields, while the action of large-scale structures of these non-linear condensations is due to a shift of 
the local background density. By definition, the Lagrangian distribution of nascent haloes of mass 
M and formation time Tf is given by 

^(q|M,r f ) = lim 6 E (q,r|M,T f )^(q,r) = &fr(q|M, r f ),5 gi (q), (21) 

where &Q(q|M, Tf) is the Lagrangian halo bias and ^'(q) is the mass density fluctuation linearly 
extrapolated to the present time To (where a(ro) = 1) and filtered on the background scale Ro; in 
Fourier space, ^(k) = (/ g i(k, ro)// g i(k, Ti n )a(ri n )) ^-^(k, r m ). One has to recall that the second 
equality in the latter equation does not mean that 5jf (q) is proportional to <5{f(q). Indeed, the 
Lagrangian bias is in general a functional of the background density field. To understand the above 
equation, one has to recall that, at sufficiently early times, the expression for the Eulerian bias 
field obtained in linear perturbation theory becomes exact (as the linear theory gets more and more 
accurate) and, going to Fourier space, 

<F(k,r) ^ § { g (k,r) = M^ a (r)#(k) (22) 

/gi(K,r ) 

when t — Ti n ^ 1- Notice that the perturbative expansion of the matter density contrast is valid 
at sufficiently early times and/or large scales, while the validity of the expansion of the halo density 

3 To solve for the integral appearing in Eq. (14) one can use, e.g., the results contained in the Appendix of 
Ref. [33]. One can also readily check that the same expression is obtained, e.g.. in the Poisson gauge. 
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contrast in powers of 5f is based on assuming a suitably large smoothing radius Rq . Assuming these 
assumptions to hold and expanding 5f*(q) in powers of <5[f(q) 

«f«l) = £&^ (#«■))', (23) 
from Eq. (8) we deduce that at first-order and in at space 

<5f (1) (k,r) ee 6f(k,r)^ i (i)(k,T)^(l + 6L(k,r))^ i ( 1 )(k,r), (24) 

To derive Eq. (24) we have used the fact that (^'^(q) can be neglected in Eq. (8). We see that the 
use of a gauge-invariant mass density contrast valid at all scales introduces a scale dependence in the 
Eulerian bias bf due to the relativistic effects which is not present if gauge-dependent density con- 
trasts are used. Of course in the Newtonian limit kr » 1 one recovers the standard scale- independent 
Eulerian bias prediction from the Lagrangian approach, bf ~ (1 + b\) = (1 + b^ 1 (M, Tf)/a(r)). 
At second order one first write the matter density contrast in Fourier space as 

V (2) (k,r) = J /Cf i (k 1 ,k 2 ;r) ( 5g i W(k 1 ,r)^ i ( 1 )(k 2 ,r) ( 5 D (k 1 +k 2 -k). (26) 

With this position one finds 

^f {2) (k,r) = |^^^6f(k 1 ,k 2 ;r)^( 1 )(k 1 ,r)^W(k 2> r)fc(k 1 +k 2 -k) > (27) 

E &i(ki,r) b\(k 2 ,r) 1 L 
6 2 (ki,k 2 ;r) = + + -b 2 (k u k 2 , r) 

1 h L ( , , /gi(k 2 ,To) (ki-k 2 ) 1 L /gi(ki,r ) (ki-k 2 ) 



2 iV 15 7gi(k 2 ,r)a(T)/ gi (k 1 ,r )^ 2 ; /gi(ki, r)a(r) / gi (k 2 , r )fc 2 2 

+ /Q(ki,k 2 ;r), 

(28) 

fefe-r) - yk.ro)/ (k 2 ,r o) 6^ 2( Mr f ) ^ 
/gi(ki,r) / g i(k 2 ,r) a 2 (r) 

Again, the Eulerian bias &f gets a scale-dependence due to relativistic effects which arises once a 
gauge-invariant definition of the density contrasts are adopted. Furthermore, in the second line of 
Eq. (28) there appears non-local terms proportional to the first-order bias parameter; they arise 
from the first-order density field evaluated at the Lagrangian point q when expressed in terms of 
the Eulerian position x [25]. The kernel ICf is conveniently computed in the Poisson gauge (rather 
than in the synchronous gauge). In the presence of large local non-Gaussianities, we can introduce 
a gravitational potential which, at some initial epoch and deep in matter domination, reads 

0in = + /&(*£ )2 - <^ n )2 » , (30) 

with the dimensionless non-linearity parameter /J^ setting the level of quadratic local NG. In the 
case of large local non-Gaussianities, /^l ^ 1> one finds [34, 35] 
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2 252 
2 



5 (VV) 2 + W 3 V,ij + 7<p' l V 2 y>,i - /^ L C -VV, 



floe 



where V 2 ^ = — ^[(X7 2 ip) 2 — <p,ik ( P'' tk ]- Inserting these expressions into Eq. (13), we obtain 



(31) 



/Cf^k^r) 



5 | 2 (ki-k 2 ) 
7 7 ^\ ^2 



18* 1 



7 k 2 



[(ki • n) + (k 2 • n)] 



+ 6i 



(ki • n) (k 2 • n 



k 2 T 



+ 



k 2 2 r 



T 



_8_ 2 
21 + 7 



(ki • k 2 

k 2 k 2 
"-l"- 2 



+ 6/; 



(ki • k 2 )' 

t 2 ^ 2 
"a ^ 2 

loc_2_ 
NL i2i.2 2 
1 2 



1 o. f loc[(ki-n) + (k 2 -n)] 



(32) 

1,2). Notice that 



PPt-3 ' 

where fc = |ki + k 2 | and we have performed an expansion in (fc^r) -1 <C 1 (i 
in the kernel the primordial non-Gaussian piece coming from the second-order density contrast is 
post-Newtonian and is damped by two powers of (hr) with respect to the Newtonian leading terms. 
The Newtonian part of the kernel does not coincide with the one of the matter density contrast 
in the Poisson gauge which can be found in Ref. [34]. Indeed, it gets a correction arising from 
the last term of the expression (13). Furthermore, there are terms which are damped by only one 
power of (fcjr); they originate from the velocity contributions in gauge-invariant definition of the 
matter density contrast and they are absent if a gauge-dependent definition of the matter density 
contrast is adopted. The same holds for the last contribution, damped by three powers of (for). 
This term comes from the primordial NG term in i/ 2 ) (which gives the dominant contribution to 
the second-order redshift perturbation (15)). All other relativistic effects have been neglected. 

The Lagrangian bias factors ^(rlM, Tf) and 6q 2 (t|M, T{) are those computed in Refs. [32] 
(see also Eqs. (16) and (17) of Ref. [25]) through the extended Press-Schechter approach and the 
peak-background split method (of course one can use refinements of the Press-Schechter halo mass 
function like the one in Ref. [36]). There are two key points in the Press-Schechter theory. On one 
side, the comoving number density of collapsed haloes is computed from the statistical properties 
of the linear density field, assumed to be Gaussian. On the other side, a patch of fluid is part 
of a collapsed region of radius R if the smoothed linear density contrast on that scale exceeds a 
suitable threshold value 6f computed at the formation time Tf. In our gauge- invariant formulation 
the first point remains of course true; the second point requires though some comments. Indeed, 
the threshold value for the matter density contrast is computed according to the spherical collapse 
model. As we have commented, the latter requires to work in the comoving-orthogonal gauge and 
therefore 6f = 5 c (t{) ~ 1.68/a(rf) corresponds to the threshold matter density contrast in that 
gauge. This does not coincide with the gauge- invariant threshold density contrast. However, it is 
easy to convince oneself that the gauge-invariant threshold value obtainable through the relation 
(10) differs from the one in the comoving-orthogonal gauge by a factor 10 -4 even for very massive 
haloes, M ~ 1O 15 M (see also the discussion in Ref. [37]). 

Adopting a gauge-invariant expression for the density contrast brings two new and interesting 
features in Eq. (28) for the Eulerian bias bf ■ First a new scale dependence parametrized by / g i(kj, r), 
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accompanied by a characteristic kernel Kf. Second, the primordial NG introduces a dependence 
on the line of sight n which comes from terms like (v • n) , which are necessary to realize the gauge- 
invariant definition of the matter density contrast. Therefore in the Eulerian bias there appears 
scale-dependent contributions which get also modulated. Notice that such v • n contributions do not 
represent the usual redshift-space distortion effects, but rather they appear as effective corrections 
to the Eulerian bias Obviously, when computing the power spectrum of two objects with the same 
bias on large scales, this modulation disappears being the power specturm real. Nevertheless, this 
does not happen when computing the power spectrum of two objects with different bias. In such 
a case, the resulting bias in the presence of some primordial NG is not only scale-dependent [22], 
but also depending on the angles cos 9 = k • n between the vector k and the vector indicating the 
line-of-sight. A computation similar to the one leading the a scale-dependent bias when some large 
local NG is included leads to a correction to the bias on large scales for two objects with different 
bias bf a (k) and bf h (k) 

Ab, (k, R U R 2 ) - WMz) k2T{k) {1 + z)k \b lA (k, R,) - b lh (k, R 2 ) j cos 6 , 

(33) 

where we have generalized our computation to a ACDM model, f2 m is the dark matter critical 
density, T(k) is the linear transfer function, /(fi m ) ~ fi^; 6 and cr 2 (R) and <J%(R) are the variance of 
the density contrast and of the velocity at a radius R, respectively. 

5 Conclusions 

In this paper we have described the computation of the Eulerian bias at second-order in perturbation 
theory. Paying attention to the gauge-invariant issues which necessarily arise when dealing with 
relativistic effects on large scales and with real observables, we have shown that some interesting 
effects show up. First of all, the Eulerian bias acquires a scale-dependence on large scales even if the 
primordial NG is totally negligible. Secondly, the primordial non-Gaussianity induces in the bias a 
modulation with the line of observation when sources with different biases are observed. Of course, 
we are well aware of the fact that our results are not complete in the sense that not all effects have 
been included. In this paper we have restricted ourselves to that part of the observed galaxy density 
contrast which is directly proportional to the dark matter density contrast through the bias. On 
the contrary, we have not discussed, for instance, the redshift-space distortion and the magnification 
effects. We leave it for future work. 
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